Coulomb corrections and thermo-conductivity of a dense plasma 
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We point out that confusion sometimes arises when using a chemical potential in plasma with 
Coulomb interactions. The results of our consideration are applied to the discussion of nuclear 
reactions screening. Finally, we present a transparent derivation of the thermal conductivity 
coefficient of a degenerate electron gas. 
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I. INTRODUCTION 

This paper considers some questions that are 
important when nuclear flames in a dense plasma are 
considered. These questions can be applied to many 
systems, for example, to supernovae of type la. If 
one considers papers dealing with nuclear statistical 
equilibrium (NSE) or with nuclear reactions screening, 
one can note that authors use different definitions of 
chemical potential and different expressions for Coulomb 
corrections. Some of these definitions are general and 
some are accurate only in particular cases. Inaccurate 
definitions can sometimes lead to the wrong conclusions 
and misleading understanding of a physical situation. 
The aim of this paper is to arrange all definitions and 
basic equations into a physical perspective and to discuss 
some paradoxes arising from inaccurate definitions. The 
second part of the article deals with kinetic properties 
of matter in white dwarfs. The question was considered 
in the 1960-1970s, [l|, 0| and references therein. One 
of the most important parameter for supernovae la 
is thermo-conductivity of degenerate electron gas. A 
method for calculating thermal conductivity coefficient 
was described in [l| but final expression was not obtained; 
whereas in a correct final expression was obtained 
but derivation is not so rigorous. There are some 
improvements in • Here we present a straightforward 
and more transparent derivation, based on ideas from [l| . 



II. CHEMICAL POTENTIAL 

A general definition of a chemical potential fj, is 
contained in the differential of internal energy E, 



dE = TdS - PdV + fidN. 



(1) 



We use standard notations here and fidN may be a 
sum for several components: fidN — > f^i^N. The 

i 

differential of the Helmholtz free energy F = E — TS 
is written as follows 

dF = -SdT - PdV + (j,dN, (2) 
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and thus it means that 

M = 



dF \ 

m) 



(3) 



T,V 



Let us consider the Coulomb free energy of a set of Ni 
nuclides: 



rik)- 



(4) 



Here is the free energy per nucleus i and it may depend 
not only on the concentration, 



Ni 
V 



(5) 



of the i-th nucleus but also on concentrations of some 
other nuclei rife. We have by definition 



Mi = 



dF n 



ON 
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(G) 



T,V k 

Then we find, for example, the pressure of nuclide gas: 
BF 

u± nuc 

T,{N k } 



Pn 



dV 



^ Mfc"fc - Fnuc/V, (7) 



which is consistent with the definition of the grand 
thermodynamic potential, f2 = — PV. 

A convenient parameter of plasma non-ideality is 
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(8) 



According to accurate Monte-Carlo (MC) simulations, 
a very good theoretical model for Coulomb corrections 
in a neutral one-component plasma is Wigner-Seitz 
approximation (in the first order of T), see references, 
e.g. in [6|. In this approximation, the correction to the 
free energy of a system is 

AF C = -0.9NiTikBT, (9) 

where a is an average distance between ions of type i, 



Ail 

4tt Ni 



1/3 



(10) 



Such a definition of at is often used (see e.g. Q), but we 
point out in advance of full discussion that this is correct 
only for one-component plasma. 
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A good correspondence of this simple model with MC 
simulation one may find, e.g. in the results by [8], where 
the main correction to free energy is approximated by 
/ = aT + ... and a = -0.899172. 

So in the case of Wigner-Seitz approxomation the 
Coulomb correction to chemical potential from equations 
© and © is 



Afi c = -1.2 



Zfe 2 



= -1.2Tik B T, 



(11) 



the result which has coefficient that is 4/3 times larger 
than the Wigner-Seitz's 0.9 because a; depends on JVj 
too. Here and below we will focus our attention on the 
numerical coefficient in the correction; thus the coefficient 
is very important. 

Sometimes <Zj is introduced in another way, through 
electron density 



a, 



3_Zi 
An n„ 



1/3 



(12) 



and the resulting expression for the Coulomb corrections 
is 



AF 



c 



-0.9NiZf /3 e 2 | ^-n 



An 



1/3 



(13) 



At the first glance such an expression leads to the 
coefficient 0.9 in (|11[) . but in this case n e depends on 
the number of ions. And if we take out one ion from our 
system then, for preserving electro-neutrality, we should 
take out Z electrons; hence, n e changes so that we obtain 
the correct answer (fTTj) . 

These considerations have direct implications for NSE 
which is essentially reduced to a Saha-like equations 
for abundances of nuclides. A correct form of these 
equations should be derived from correct expressions for 
chemical potentials in equilibrium. However, it is also 
used sometimes in another approach, when Coulomb 
corrections (with the coefficient 0.9) are introduced 
directly to the binding energies of reacting nuclei. This 
approach also seems natural, at least while nuclides form 
an almost ideal gas (see, e.g., a widely cited paper [9j), 
but the results for nuclear abundances will differ from 
those obtained with the correction (fTTj) for chemical 
potential. 

We can derive the result (fTTj) from a more general 
expression. If Qi denotes the Coulomb correction to the 
free energy and has the functional form 



Qi = k B Tf (Ti 



(14) 



then we find the Coulomb contribution to the chemical 
potential as 



Coul 



fo(Ti 



1 df 



3<91nr,; 



(15) 



this agrees with (fTTj) in case of the Wigner-Seitz 
approximation. But we want to point out that it 



contradicts equation (17) of Q often used in NSE studies. 
Our expression (|15[) does not coincide also with equation 
(97) of the famous review by [lfj- It was fount in (lo| 
that fii — fcsT/o(rj), but this expression leads to zero 
correction to P nuc in equation ([7]). On the other hand, 
our equation (|15p gives the expected correction to P n uc 
after substitution into equation (f7j). 

It is a very important point in our reasoning, so let us 
consider it carefully. As will be shown in section IIII[ a 
more appropriate expression for Coulomb corrections is 
(fTJJj), and not ©-©■ But if we use the assumption of 
electro-neutrality, the numbers of electrons and ions are 
directly related. It is correct in case when we consider our 
system on a timescale that is much larger than plasma 
frequency (like in NSE) - so we cannot separate electrons 
and ions. In this case, the chemical potential is defined 
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(16) 



electro — ncutr. 



and for one-component mixture, the answer is (|11[) . Here 
the full correction could be ascribed to ions if electrons 
are considered as background or to electrons if ions are 
background. But when we consider kinetic properties of 
our medium, we study it on timescales when electro- 
neutrality could be violated — so we should consider 
electrons and ions separately, like the corrections 
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these corrections should be applied to (]4l)j) below (there, 
in a kinetic case, all thermodynamic quantities have an 
approximate meaning, and have to be carefully defined 
similar to that in [11(). Here we should emphasize one 
important point. Coulomb corrections in dense plasma 
arise due to electron-ion interactions because all ions are 
screened by electron clouds. The notation Fa that is used 
in many works is a bit misleading because it is not pure 
ion-ion interaction. 

In the book by [l2| the authors rely mostly on free 
energy F instead of fii. But in many cases correct 
knowledge of /x is required. Actually, the expression which 
is equivalent to (|15p appeared already in [13]. In the paper 
[l4j . the result for the \i coefficient is close to 1.2, but it 
was obtained using Monte-Carlo sumilation and rough 
expansion. However, the argument on the role of the 
chemical potentials and debates on approaches to them 
still persists. For example, in [ll| the difference of their 
use of chemical potentials have been discussed with that 
of 0. Essentially, in 1 15] it was found that the Coulomb 
correction found in [l4| for a classical plasma is applicable 
to quantum plasma as well. (See also [lj| for a detailed 
review on modern perspective on Coulomb systems.) 

Since there is some controversy on the form of Coulomb 
corrections in the literature, in the paper (l7| the results 
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for the Coulomb correction from equation (fT4")l . as was 
taken by |9(, with the effect of the additional term 
in equation (fT5|) have been compared. The effect of 
the new term in (|15[) on the proton fraction (which 
is very important in collapse computations) is visible 
only at highest density but still very small. A bit more 
pronounced effect is obvious for the average mass of 
heavy nuclei again at highest density when the NSE 
model is near the border of applicability. However, a more 
detailed studv_is needed here along the lines undertaken 
recently by 
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III. SCREENING OF NUCLEAR REACTIONS 

A very important application of the previous reasoning 
is in the field of a classical nuclear burning screening. 
Consider the regime of strong screening T > 1 . According 
to [HI the expression for screening factor in this case can 
be written as 

cxp (#i 2 (0)) = exp(- {^Zt + Z 2 ) - fxiZt) - fi(Z 2 )} /T). 

(18) 

Here and below ks = 1- So theoretically, using the above 
definition of a chemical potential, the screening factor is 



#12(0) = 1.2 (r Zl _ 



(19) 



In order not to create confusion, let us consider a 
commonly examined situation when Z\ = Z 2 = Z . We 
obtain in the leading order (here we will discuss only the 
leading order of corrections) 



#i 2 (o) = i.4ior z , 



(20) 



(here 1.410 = 1.2(2 5 / 3 - 2)). On the other hand if 
we consider the definition of screening factor using free 
energy, 



H 12 (0) = -(F(2Z)-2F(Z))/T, 



and using ^ obtain 



#i 2 (0) = 0.9 {T 2Z - 2T Z ) = 1.05^ 



(21) 



(22) 



In the case of T 3> 1 the difference in (fT9")l and (|2"2")l can 
be extremely high. In some works the chemical potential 
has coefficient 0.9 and so using two definitions of #12 the 
answer (|2"2")l is obtained. 

First of all, let us consider the physical sense of a 
chemical potential. Equation (|18[) tells us that we take 
out from our system two ions with a charge Z and put 
back one ion with a charge 2Z. From the chemical point 
of view, ions with Z and 2Z are different elements and the 
average distance between Z-ions changes. So if we recall 
that coefficient 1.2 arises from the change of ai, then (|2"2")l 
is incorrect. But for finding a proper screening correction 
one should stand on electromagnetic force point of view: 
these are elements of the same nature - and there should 
be only one distance between reactants. 



At this point a reader can object in the following way: 
a nuclear reaction is a fast process, but the effect we are 
discussing now exists because of a global (on scales of our 
system) variation of the average distance between ions: 
this is a slow process. Can it really change the repulsion 
potential — a barrier, our ion is penetrating through? 
Theory of such corrections rely on a statistical averaging 
[1J, [19J and operate with average quantities, so in such a 
theory we should consider these effects. 

On the other hand, if we consider definition of free 
energy using n e (|13p . in nuclear reactions, the electron 
number density n e stays constant, unlike the situation 
when we define chemical potential and take one ion with 
its electron cloud from our system. The average distance 
a,i does not change in these terms, and then the coefficient 
0.9 should stay in (|19p . We see that there is some paradox 
in this situation: one reasoning leads to 0.9, another to 
1.2. 

Let us resolve the paradox. Expression (f2"2l used 
by scientific community is correct. The reason of 
the misunderstanding is the definition of which is 
sometimes incorrect. What is the Wigner-Seitz model? 
This model tells us that the most significant Coulomb 
corrections are due to the interaction of ions with 
electrons (and defining these correction through the 
average distance between ions may result in a confusion 
of the reader). Electron gas is uniformly spread over 
the whole volume (we omit such effects as electron 
polarization etc) and we should define aj for an ion Zi as 
a radius around the ion, which contains electron charge 
(-Zi): 



47T ?l e 



1/3 



(23) 



Here, in general, is not the average distance between 
ions: this would be true only if we had a one-component 
mixture - what is undoubtedly wrong in nuclear 
reactions, our mixture after reaction is multicomponent. 
So the correct definition of at is a radius of electron cloud 
with charge (—Zi). Such an electron coat screens nuclear 
repulsion potential and one ion from another — there is 
no direct Coulomb interaction between ions. As we have 
mentioned above, Ti e is constant, so the correct leading 
order in screening is 



#12(0) = 1.057T* 



(24) 



Finally we could write the full classical correction from 
MC simulation with equal charges (see e.g. [8j) 



1.039957r^ 323064 



#12(0) = 1.056299r z 

-0.545823 lnr z - 1.036319. (25) 



Let us emphasize one moment. When defining chemical 
potential, we should obey the electro-neutrality law and 
take one ion from the system with its electron cloud. And 
in nuclear reactions nothing is taken from the system, 
so the chemical potential has another physical meaning, 
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and the correct definition is done by free energy (PHI) . 
A very bright demonstration of this point is the famous 
review by Ichimaru @ with its section C on "Chemical 
potentials": only free energies are discussed, and fi symbol 
does not appear there. 



IV. THERMAL CONDUCTIVITY 

It is commonly supposed that write dwatfs (WD) are 
progenitors for supernovae la. On the initial stage of 
explosion, thermo-conductivity plays a significant role as 
it drives the flame front of thermonuclear burning [20l |. 
Typical parameters of matter in presupernovae WD are 
T ~ 10 9 K and p ~ 10 9 g/cm 3 . Chemical composition 
is a mixture 12 C+ 16 0. For these conditions the electron 
gas is relativistic and degenerate ep 3> fcsT, pf 3> m e c, 
and the parameter of non-ideality © T ~ 1. 

Let us consider an electron gas which is described 
by its distribution function in a phase space /(p,x, t). 
Boltzmann equation for evolution of such a function is 
the following (in this section we will use units in which 
c = h = ks = 1 and in the end we will restore 
dimensionality of variables): 

| = St/ = - J W (p, pO [/(I-/')- /'(I -f)]^. 

(26) 

Here we assume that the scattering of electrons occurs 
only on ions, and the distribution function of the latter 
does not change. Under our conditions the ions arc 
non-degenerate, non-relativistic and have large mass in 
comparison with electrons mass; they can be treated as 
being at rest. An approach which considers the scattering 
on one ion (separately from other ions) is incorrect when 
r ^> 1 because of the correlation in ion positions. Thus 
we shall consider a situation when T < 1 and suppose 
that our answer has a right limit r — > 1 (finally we want 
the answer for T ~ 1). This is true because when T 1 
we have an ideal gas behavior of ions. The case T ~ 1 is 
called Coulomb liquid: there is no significant correlation 
in positions and full correlation sets in the solid state 
when T - 170 (see [UHl]!). 

Cross-section of scattering on an ion at rest (with 
potential Aq{t) = Ze/r) is written as follows [24| (we 
will use a short definition dp' = 2d 3 p' / (2ir) 3 ^" 

da_ f " 2 

dp 1 

where Ao(q) = 47rZe/q 2 . The S- function tells us that the 
scattering is elastic. The comparison of e-i and e-e cross- 
sections gives a e i/a ee ~ Z 2 3> 1, see e.g. [24j |. so we omit 
the e-e scattering here. 

As usual, we consider a small deviation from the 
equilibrium distribution function: 



Z. = ^5{e p , - e p )\A (p' - P)| 2 (l - |5r) g|, (27) 



1 



f P = fp-^P = / P O +^/ P °(l-/p O )0 P , /> 



T 



1 



So Boltzmann equation in the case of thermo- 
conductivity looks like 



df° e — w 



(W)T=-i J w(p ,p')(^0')/°(l-/ )^. 



(29) 

In the case of degenerate electron gas, sT <C /i, so 
we could substitute /i instead of in. Let us define 
the abbreviation T p ^ p , = w(p, p')/°(l - /")• After 
multiplication on <fi and integration over the whole p, we 
can reduce the Boltzmann equation to the form 



7 2 

JE 



1 



kT 2T 
where energy flux is defined as 



2 2d 3 p 2d 3 p> 



(2tt) 3 (2tt) 3 : 



Je = |(/-/>(e-M)-0, 



(31) 



and the thermo-conductivity k is defined in the following 
way: 



JE 



(32) 



Equation ()30|) has exact physical meaning - it can 
be shown (see |25j) that the rhs of (|3U|) is a rate of 
entropy generation by particles' scattering. On the lhs 
is the entropy generation rate through conductivity. The 
following method of coefficient determination is good - 
we could substitute arbitrary functions (f> and get answers 
for k. There exists a theorem that the maximum value of 
k for different <f> is mathematically exact. 

This is a simple case, namely, isotropic. There is a 
limited set of quantities in our task; that is why we could 
guess the answer 



A(e p - fi)p a d a T. 



(33) 



The previous expression can be simply explained: it 
should be proportional to the temperature gradient, 
the only physical quantity that can create a scalar 
with combination with gradient is a momentum. The 
factor (e p — fi) was introduced because all macroscopic 
parameters should be defined through /q but not through 
the correction proportional to (p. The factor A is 
a dimensional constant, but it will be cancelled out 
immediately after the insertion of <f> into (|30p , so we omit 
it now. 

Let us consider the integration 



J'e = ^ |(e p -/i)v/ p °(l-/ p >dp. 



(34) 



Choosing for definiteness d a T = (0, 0, VT) and using 
the fact of high degeneracy of electrons (therefore 



I m d ~tde = - f TW) ~ S^V'V) + •••) 



o 

we obtain 



(e p -Ai)/T _|_ i ■ 

(28) 



^ = g-(VT)TV-m 2 ) 3 / 2 . 



(35) 
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The rhs of equation ([50)) : 

mis : / (d> - 0') 2 r 



after all substitutions 



2d 3 p 2d 3 p' 
" v ' (2tt) 3 (2tt) 3 ! 



(36) 



RHS = /(e p - M) 2 fc(^r)„(^T) TO 2 5( e - e')|A | 2 (l - ^) ^ • -/ p (l - f°)mv ^ g*V . 
Doing the same integration, 

RHS = ^e^T 3 ^ 2 - m^/^VT) 2 ! 1 9 3 |A ( g )| 2 (l - dg. 



(37) 



(38) 



r 



Let us denote Coulomb logarithm (recalling that 
A (q) = AnZe/q 2 = 4nZeip{q)) 



2p 



Ki = / g 3 1 </?(<?) I 



(39) 



RHS = ^Z 2 eVTW - m 2 ) 1 / 2 (VT) 2 l A el . (40) 
9 p 

After the substitution of ([35]) and (gOJ) in ([301), and 
defining vf = p/e, we obtain 



1 T(/i 2 -mg) 5 /y 
367r Z 2 e 4 ni/iA e i ' 

recovering dimension 

1 fc|T( M 2 - m 2 c 4 ) 5 /V 



(41) 



(42) 



367r (he) 3 Z 2 e 4 ni)iA e i 

This is the same result as calculated in Q. 

Last uncomputed quantity is the Coulomb logarithm 
Because of quantum nature of our scattering (see 



leg I 
fit) 



Ze 2 



= « 1, 



(43) 



then the minimal momentum transferred is g m j n = H/l, 
where I is the typical length of interaction. Then (13TH) can 
be rewritten as 



L, 



In 



2 Pe l 



2c 2 ' 



(44) 



We should define I - maximum length with interaction. 
Recalling reasoning from the previous paragraph we 
could put I = a,i, so in this case the Coulomb logarithm: 



In ((IStt) 1 /^ 1 / 3 ) - || = In (3.838Z 1 / 3 ) 



V F 

2c 2 ' 

(45) 



In hydrodynamics there are general expressions for mass 
and heat fluxes (see, e.g. |26]): 



i = -aVfi-/3VT, 
q = -<5V/i-7VT. 



(46) 



Here we use a definition of q with subtraction of Let 
us restore terms with V/i in kinetic equation: 



df e ~ n df , s 

_ — (vV)T+— (vV) M = 



w ( P ,p')(0-0')?rrfp'- 

oe 



Like above we look for solution in the form 

<j) = A{e- M)(pV)T + £(pV)M- 

And using the same techniques, we obtain 

4> = A(e-//)(pV)T + AT(pV)/i, 

A = 1 (/i 2 - m 2 e )v F 

4tt Z 2 e 4 "n i fxA ei T' 

Based on definitions 



(47) 
(48) 

(49) 



q = - j (e - n)v^cpdp' , 



(50) 



we obtain Onsager's result S = T(3 and in our case (of 

2 

degenerate electrons) a = —AT^j-^,] j3 = 0; 8 — 0; 7 = k. 

At the end we can give some limits of applicability for 
our expression for k. As we can recall, we have used 

00 

m^de = ^M-^T 2 f( /1 )-^TV"'M+. • • 

(51) 

in our calculations. We have omitted the 3rd term and 
this is correct when 



T < 



' 60 ?/>"(M) 
7tt 2 ^""G") 



0.16/i, 



(52) 
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where the last equality arises from our previous 
calculations. This is the validity criterion for our 
expression (|42j) . 

V. CONCLUSIONS 

This work presents methodical notes on Coulomb 
corrections in a non-ideal plasma. A quantity cti that 
enters the corrections in different works is defined in a 
different manner, what can lead sometimes to inaccurate 
results. We emphasize that a correct definition of at is a 
radius of electron cloud with charge (— Zi): 

— (s^) 1 "- <53 » 

This quantity is used for definition of free energy 
When the relation (1531 is used the physical sense of 
correction to F is clear: they arise due to interaction of 
ions with electron cloud, and not ion-ion interactions. 

An example of quantity that depends significantly on 
such corrections is classical screening of nuclear reactions. 
Sometimes chemical potentials are used for definition of 
#12(0). We argue that a chemical potential does not 
reflect the essence of the matter: the interaction is fully 



electromagnetic in a screening process, and elements with 
different Z have qualitatively the same nature relative to 
these interactions, while no particles is taken from the 
system in nuclear reactions. A correct formula should be 
written in terms of free energy and should read 

#i 2 (0) = - [F(N - 2, f ) - 2F(N, 0)] /T, (54) 

where F(N — 2,1) is a free energy for N — 2 ions with 
charge Z and one ion with charge 2Z, while F(N, 0) is 
defined respectively for ions before reaction [l9j . 

Finally, we give a straightforward derivation for a 
thermo-conductivity coefficient of degenerate relativistic 
electron gas. 
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